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Abstract 

Recently, a system with uniformly hyperbolic attractor of Smale- Williams type has 
been suggested [Kuznetsov, Phys. Rev. Lett., 95, 144101, 2005]. This system consists 
of two coupled non-autonomous van der Pol oscillators and admits simple physi- 
cal realization. In present paper we introduce amplitude equations for this system 
and prove that the attractor of the system of amplitude equations is also uniformly 
hyperbolic. Also we represent qualitative illustrations as well as quantitative char- 
acteristics of a chaotic dynamics on this attractor. 
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1 Introduction 



One of the interesting problems in nonlinear dynamics is a development of 
physical applications of the hyperbolic theory. A trajectory in phase space of 
a dynamical system is hyperbolic when a tangent vector space on each its point 
can be described as a direct sum of subspaces spanned on stable and unsta- 
ble perturbation vectors and this representation is invariant along the trajec- 
tory. Dissipative system with attractor containing only hyperbolic trajectories 
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demonstrates strong chaotic properties and permits advanced mathematical 
analysis. Monographs and textbooks on nonlinear dynamics provide examples 
of such attractors, but most of them are artificial mathematical constructions 
like Plykin attractor and Smale- Williams solenoid [2^|B|7f2TpiJl[2lfT7ll5] . On the 
other hand, attractors of realistic systems with complicated dynamics, like the 
Lorenz model [X1T16] . do not relate to the class of uniformly hyperbolic, so the 
hyperbolic theory can not be applied to them in corpore. 

In a recent paper [13] one of the authors suggested an implementation of the 
Smale- Williams attractor in a system of two coupled non-autonomous van der 
Pol oscillators with natural frequencies relating as 1:2. Subsystems are acti- 
vated alternately due to a slow variation of their excitation parameters by an 
external force. The excitation is passed from one oscillator to another, so that 
a phase of the oscillations is doubled within each full cycle. This system was 
realized as an electronic devise and studied experimentally |12j . Assumption 
about the hyperbolicity of the attractor was based on the fact that a phase 
of one of the oscillators, being measured at successive stages of its excitation, 
was obeyed to the Bernoulli map, as must be for the Smale- Williams attractor. 
Moreover, computations revealed a robustness of the attractor: Its Cantor-like 
transverse structure and positive Lyapunov exponent were insensitive to vari- 
ations of parameters in the equations. Further analysis [11] showed that 4D 
phase space of a Poincare map, describing a successive states of the system 
over a period of the external forcing, contains a toroidal absorbing domain (a 
direct product of ID circle and 3D ball) that is mapped into itself according to 
the Smale- Williams procedure. Performed numerical analysis proved the va- 
lidity of sufficient conditions of the hyperbolicity that are formulated in terms 
of expanding and contracting cones in tangent vector space. 

An idea of alternate excitation of two oscillators that pass the excitation from 
one to another can be applied to amplitude equations which are derived for 
the van der Pol equations via the method of slow varying amplitudes. Actu- 
ally, equations of this type can be obtained for a wide variety of dynamical 
systems because, in fact, they correspond to a normal form of Andronov-Hopf 
bifurcation of birth of a limit cycle. For example, amplitude equations are 
employed in Landau's theory of turbulence [ToTTi] . 

Reasoning in the manner of Landau's theory, we can assume that a spatially 
extended system have two modes with frequencies relating as 1:2. Some addi- 
tional periodic and slow component acts on these modes so that they excite 
alternately and pass the excitation from one to another through a nonlinear 
coupling. Then, according to the idea of [13] , this can result in the formation 
of a hyperbolic attractor of Smale- Williams type. 

This paper is devoted to a study of a system of two coupled non-autonomous 
amplitude equations. We provide a proof of hyperbolicity of an attractor of 
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this system and also discuss some attributes of its hyperbolic dynamics. 



2 System of two coupled non-autonomous van der Pol oscillators 
and corresponding amplitude equations 



The starting point of our analysis is a system of two coupled van der Pol 
oscillators that has been proposed in [13J: 



c 2 ]x + UqQ 

y - [-A cos(27rr/P) - y 2 ]y + Au 2 y = ex 2 



x — [A cos(27rr/ P) — x }x + u x = ey cos ujqt, 

,2„. _ _2 



where r denotes time. The oscillators have natural frequencies io§ and 2u; , 
respectively. Their bifurcation parameters controlling a birth of a limit cycle 
undergo a variation with a period P and amplitude A. The variations are in 
counter phases: when the first oscillator is excited the second one is not and 
vice versa. The forcing is supposed to be slow, i.e., its period is much larger 
than natural periods of the oscillators. The oscillators are coupled and the 
intensity of this interaction is controlled by a parameter e. 

The subject of our study shall be a system of amplitude equations correspond- 
ing to Eqs ([1]). To derive these equation we follow the standard technic and 
assume a solution to Eqs ([T]) to be oscillations with frequencies u and 2u , 
respectively, and with slow varying complex amplitudes: 

x(t) = a{r)e iuJ0T + a*(r) e -^° T , 

y ( T ) = b(T)e 2luJ0T + b*(r)e~ 2iuJ0T . ^ 

Upper index "*" means complex conjugation. The derivatives of complex am- 
plitudes a and b should satisfy additional conditions 



ae lu30T + a*e- lu)0T = 0, 

be 2iuJ0T + b* e - 2iuJoT = 0. 



(3) 



Taking these conditions into account we substitute (j2J) into §I§ and average 
over the period 27i/u . Because the complex amplitudes a and b are supposed 
to be slow, the resulting equations have the form: 

d = (A/2)acos(27rr/P) - \a\ 2 a/2 - ieb/{4u ), 
b = -(A/2)6cos(27rr/P) - \b\ 2 b/2 - tea 2 /(4uj ). ' ' ' 

These equations allows the following rescaling of time variable and parameters: 

t = r/2, T = P/2, e = e/2to . (5) 
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As a result, we obtain a sought system of amplitude equations: 



b 



a = Aacos(2irt/T) — \a\ 2 a — ieb, 
• = -Abcos(2nt/T) - \b\ 2 b - iea 2 . 



(6) 



To clarify the nature of dynamics of the system ([6]), we consider a behavior 
of amplitudes and phases of complex variables a and b. Define phases within 
the interval [0, 27r): = arga, ip = argfr. Suppose that the first oscillator is 
excited and its amplitude \a\ is high. Hence, the second one is suppressed so 
its amplitude |6| is small. The coefficients in Eqs. ([6]) are real except the cou- 
pling. It means that the phases can be changed only as a result of interaction 
between subsystems. But when a is excited |6| is small and its action on a is 
negligible. Thus, the phase of a remains constant during the excitation stage. 
The opposite influence from a to b is high and the coupling term is propor- 
tional to a 2 . It means that after a half period T/2 at the threshold of its own 
excitation the oscillator b inherits a doubled phase of a (the phase also gets 
a shift — 7r/2 because of the imaginary unit at the coupling term). Now the 
roles of the subsystems are exchanged. The phase of b remains constant when 
this subsystem is excited and at the end, after the other half period T/2, the 
phase is returned back to a through a linear coupling term (also with the shift 
— 7r/2). As a result, the first oscillator a doubles its phase during the period T. 

Fig. [1] demonstrates a numerical solution to Eqs. (jfJJ) illustrating these consid- 
erations. Here and below the solutions are obtained via Runge-Kutta method 
of the fifth order with a step-size control [TU] . Panel (a) demonstrates an alter- 
nate excitation of the subsystems a and b. Notice, that each time the amplitude 
\a\ attains its maximum, the height of the maximum is a bit different from 
the previous one. The same is for the amplitude |6|. This is a manifestation 
of a chaotic nature of the discussed dynamics. Panel (b) represents a time de- 
pendance of a complex variable a. We observe a series of spikes corresponding 
to the stages of the excitation of a. Projections of spikes on the plane {Re a, 
Ima} are almost straight lines. This confirms our previous conclusion that the 
phase of an excited subsystem, i.e., an angle between the projection and the 
real axis Re a, remains constant or, at least, almost constant. This panel also 
illustrates the phase doubling after each period T: compare angles between 
spikes 4, 5 and 6. 

The above discussion allows us to write down a map for a series of phases 
4> n = arg a(nT) that are measured over the time step T: 



Up to a constant addition (that can be eliminated by a shift of the origin 
of the phase) this map coincides with the well known Bernoulli map [20]. It 
demonstrates a chaotic dynamics and the chaos is homogeneous: a rate of 



4> n+ i = 20, 



n 



7T mod 2tt. 



(7) 
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Figure 1. Illustration of a dynamics of the system ([6]): alternate excitation and 
decay of the components a and b (panel (a)) and time dependance of real and 
imaginary parts of a (panel (b)). Observe that being excited the component a has 
almost constant phase, but the phase is doubled at every new stage of the excitation 
(compare, for example, spikes labelled 4, 5 and 6). 
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Figure 2. Numerically computed map for the phase <p n = arg a(nT) measured over 
the time step T at A = 3, T = 5 and e = 0.05. 

exponential divergence of two close trajectories is identical in each point of 
the phase space and equal to In 2. 

Fig. [2] show 4> n+ i = arg a{nT + T) vs. <p n = arg a(nT) computed numerically 
for the system (|6]). One can see a very high correspondence with the map ([7]). 
The most important point is the topological equivalence between the empirical 
map and the Bernoulli map: one full circle that passes a preimage (i.e., turn 
of a phase on 2tt) implies two passes of a circle for an image (47r turn of a 
phase). 

The correspondence with the Bernoulli map presumes that among Lyapunov 
exponents of the system Q the one should be equal to In 2/T. To compute the 
Lyapunov exponents we employ the Benettin's algorithm [51TT8] that requires 
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to solve simultaneously Eqs. ([6]) and four sets of linearized equations for small 
perturbations: 



5xq = Acos(2iit/T)5x — (3xq + x\)5xq — 2xqX\8x\ + eSx 3 , 
5ii = Acos(2nt/T)5xi — (xq + 3x\)5xi — 2xoXi<5xo — e#x 2 , 
6x2 = —Acos(2iit/T)5x2 — (3x1 + x t)^ x 2 ~ 2x 2 x 3 5x 3 + 2e(x 5xi + Xi<5x ), 
5x3 = —Acos(27it/T)5x3 — (x\ + 3x3)5x3 — 2X2X3&E2 — 2e(xo5xo — x\5xi), 

(8) 

where xq = Re a, X\ = Ima, X2 = Re b, X3 = Imb, and 5xo, 5x\, 6x2, SX3 
denotes small perturbations to these values. Before the start of the compu- 
tations, we initialize each set of linearized equations by a unit perturbation 
vector, so that four these vectors comprise an orthogonal system. Advancing 
the solution, we perform the Gram-Schmidt orthogonalization [18J of these 
vectors after each time interval T and normalize them to prevent a numerical 
overflow. Four Lyapunov exponents Xk appear as time averaged logarithms of 
norms of the perturbation vectors. Finally, we multiply the computed values 
on T, to obtain the Lyapunov exponents for corresponding map, describing 
the dynamics of the system with time step T: = A&T 



Fig. [3] demonstrates the exponents as functions of the parameter A. The 
largest exponent is positive for a wide range of parameter values. This is the 
evidence of a chaotic nature of the observed dynamics. Three others are always 
negative. For example, at A = 3 the Lyapunov exponents are Ai ~ 0.691, 
A 2 ~ —4.06, A 3 m —6.48, and A 4 rs —9.06. As typical for maps and non- 
autonomous continuous time systems, a zero exponent is absent. The largest 
exponent is close to In 2. This corresponds to our qualitative conclusion that 
the dynamics of the phase of a is obeyed, to some degree of accuracy, by the 
Bernoulli map ([7]). Also notice that the largest exponent is almost independent 
on A, and the others vary rather smoothly, without peaks and canyons. This 
can be treated as a manifestation of a robustness of the observed chaotic 
oscillations which, in turn, indicates a hyperbolic nature of an underlying 
chaotic attr actor. 



This behavior of Lyapunov exponents agrees very well with the behavior of 
the exponents of the initial system of two van der Pol oscillators [T3][l~2"] . 
In Fig. [3] the Lyapunov exponents of the initial system are plotted by dashed 
lines. Values of parameters, being recalculated as required by Eq. (jSJ), corre- 
spond to the parameters of the system (jBJ). The largest Lyapunov exponents 
of two systems are almost identical and three others pairs are also close to 
each other. We treat this good correspondence between the initial system and 
the approximate amplitude equation as one more demonstration of a robust- 
ness. In turn, the robustness is related to a hyperbolic nature of the observed 
dynamics. 
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Figure 3. Lyapunov exponents vs. parameter A for the system (j6|) at T = 5 
and e = 0.05. Dash lines indicate corresponding Lyapunov exponents for the initial 
system ([1]) at P = 10, ujq = 2ir and e = 0.2ir. 

3 Stroboscopic Poincare map and absorbing domain 

Consider our system in a certain time moment t n = nT. Its instantaneous 
state is given by a vector x n = {Rea(t n ),Ima(t n ),Reb(t n ),lmb(t n )}. If we 
take this vector as an initial state and integrate the equations over the period 
T, we find a new vector x n+1 = {Rea(t n+ i), Ima(t n+ i), Re6(i n+ i), Imfe(t n+1 )}. 
This vector is unambiguously determined by the vector x n . Hence, we can 
define a mapping that operates in M 4 : 



This equation defines a stroboscopic Poincare map for the system ([6]). From 
a geometric point of view, each vector x n belongs to a 4D hyperplane t = 
t n = nT which is a section of a flow of trajectories in the 5D extended phase 
space {Re a, Im a, Re b, Im b, £}. Because of periodicity of the phase in t, these 
hyperplanes can be identified. So, we can say that Eq. Qj maps 4D hyperplane 
{Re a, Im a, Re b, Im b} onto itself. 

The map T is generated by differential equations whose right-hand sides are 
smooth and limited within a finite region of a and b. According to the existence, 
uniqueness, continuity and differentiability theorems, the map T is one-to-one 
diffeomorphism belonging to class C°°) pE]. 

In course of iterations, the map x n+1 = T(x n ) expands a volume element 
in a direction, associated with the phase 0, see Eq. (J7J), while three others 
directions are contracting. Taking into account a periodicity of the phase, we 
can consider a 4D toroid which is a direct product of a ID circle and 3D ball. 
One iteration of the map corresponds to a longitudinal stretch and transverse 
contraction of this object followed by its kinking and insertion into the initial 
area as a double loop. This exactly corresponds to the Smale and Williams 




(9) 
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procedure except that our map is embedded into 4D rather then 3D state 
space. 

Mentioned toroid with its interior is an absorbing domain which we denote as 
U. Points, belonging to this area, fall into the interior of U under the action 
of the map T: T(U) C Int U. An attractor may be defined as an intersection 
of all images of U, generated by successive iterations of T: A = fl^Li T n (U). 

Consider the space {Re a, Im a, Re b, Im b} and define a torus as 

{\a\-rf + \b\ 2 = { X df. (10) 

Here r denotes a radius of the ID circle, and \d is a radius of the 3D ball 
that appears in a radial section of the torus by a hyperplane perpendicular 
to the circle. Parameter x is introduced to scan points both on the surface 
(x = 1) and inside the absorbing domain, < x < 1. The torus fflUj) can also 
be defined parametrically via the equations 

Re a = (xd cos 9 + r) cos 0, Re 6 = x<isin6 l cos'?/>, 
Im a = (xd cos 9 + r) sin </>, Im b = xd sin 9 simp. 

Here parameters <fi and ip can be treated as phases of a and b, respectively. 

Numerical values of r and d should be found for the attractor to be fitted 
inside the torus. Consider specific values of the parameters of the system, say, 
A = 3, T = 5, e = 0.05, and compute a large amount of points belonging to the 
attractor of the map (Q. Projections of these points to the plane {Re a, Ima} 
fall somewhere in a vicinity of a circle whose radius r can be found as a mean 
amplitude \a\ of these points. For our specific parameters we find r = 0.631166. 

To compute d, let us consider a series of tori parameterized by the section ra- 
dius d x that grows from with some small step (we assume here that x — !)■ 
All the tori have identical radiuses r. For each d x we cover the surface of a cor- 
responding torus by a mesh, using parametric Eqs. ffTTj) . In our computations 
the angle variables <p, ip and 9 were varied with the step 27r/50. The map (Q is 
iterated once form each node of the mesh. An image point falls on the surface 
of a new torus whose section radius is equal to J(\a\ — r) 2 + |6| 2 , where \a\ 
and |6| are related here to the image point. Thus, the torus d x generates a set 
of "image" tori. Because we are seeking for an absorbing domain, we must take 
the worst case, i.e., the image torus with the largest section radius. Denoting 
this radius as d y , we obtain a function d y = F(d x ) that is plotted in Fig. HI 
This function monotonously grows and meets the line d y = d x at d x m 0.0918. 
All tori to the right from this point can be taken as an absorbing domain 
because their surfaces, as well as the interior, are mapped into their interior. 
We take a torus whose section radius d is about 1% above the intersection 
point. So, the set of parameters of the system and corresponding parameters 
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Figure 4. Selection of the parameters of an absorbing domain (see explanations in 
the text). 




Figure 5. 3D projection of a toroidal absorbing domain and enclosed attractor of 
the Poincare map Q. 



of the absorbing domain are the following: 

A = 3, T = 5, e = 0.05, 
r = 0.631166, d = 0.0927. 



Fig. [5] demonstrates 3D view of the absorbing domain with an enclosed attrac- 
tor. Points of the attractor are project from 4D space {Re a, Im a, Re b, Im b} 
along axis Imb onto the section plane ip = 0. 
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4 Sufficient condition for hyperbolicity: general formulation of a 
method 

To verify hyperbolic nature of the attractor of the system ([6]) , we shall follow 
the method developed in [TT] for the attractor of two non-autonomous van der 
Pol oscillators ([I]). In this section we reproduce the details of this method. 

The central point is the theorem on expanding and contracting cones [22][T0ll9] . 
Unlike the general case, it is sufficient here to deal with a diffeomorphism of 
class C°° in the Euclidian space M 4 , namely, we consider the Poincare map 
T(x). Let DXj. be the Jacobi matrix of the map x' = T(x) at x: DT X = 
{dx'Jdxj}, i,j = 0, 1,2,3, and let DT" 1 designates the derivative matrix for 
the inverse map T _1 (a;). Also, let 5x denotes the vector of a small perturbation 
to x. In a linear approximation the evolution of a perturbed state x + <5x 
corresponds to transformation of the vector <5x according to a linear mapping 
dx.' = DT x 5x. Vectors <5x form a tangent space Y x associated with x. 

Theorem ([22.10.9J) Suppose that a diffeomorphism T of class C°° maps 
a bounded domain U C M 4 into itself: T(U) C Int U, and A C Int U is an 
invariant set for the diffeomorphism. The set A will be uniformly hyperbolic if 
there exists a constant 7 > 1 and the following conditions hold: 

(1) The expanding and contracting cones S2 and C2 may be defined in the 
tangent space Y x at each x G A, such that \\DT x u\\ > 7||u|| for all 
u G S2, and ||DT^ 1 v|| > 7||v|| for all v G CJ; moreover, for all x G A 
they satisfy S% n Q = and S2 + C2 = Y X . 

(2) The cones S2 are invariant with respect to action o/DT, and C2 are 
invariant with respect to action o/DT _1 , i.e., for all x G A T)T X (S2) C 
Sl (x) andVT2 l (C2)cC^ 1{x) . 

If the formulated conditions are valid for all points of the absorbing domain 
U, they are obviously true for the attractor A C U. Therefore, the following 
procedure can be performed for a verification of the conditions of the theorem. 

Starting at x = {Re a, Ima, Re 6, Im6} G U, we solve Eqs. ([6]) numerically 
on the interval t G [0,T] and get the image x' = T(x). Also, we initialize 
four sets of equations for small perturbations (jSJ) with unit vectors {1, 0, 0, 0}, 
{0, 1, 0, 0}, {0, 0, 1, 0} and {0, 0, 0, 1}, respectively, and solve these equations si- 
multaneously with the original system. From the resulting four vector-columns 
we compose a matrix U = DT S . 

If the Poincare map is iterated one time from x, any perturbation vector u 
transforms to u' = Uu. A squared Euclidean norm of this vector is ||u'|| 2 = 
u T U T Uu where T means the transposition. Via the inverse matrix U _1 we can 
write u = U _1 u' and ||u|| 2 = u /T U _1 ' T U _1 u'. A condition that the preimage 
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of u' relates to the expanding cone S% is an inequality ||u'|| > 7||u||, or 

u' T H 7 u' < 0, (13) 

where H 7 = tr^lT 1 - 7 ~ 2 ; . 

If we start from x' = T(x), a vector u' transforms to u" = U'u', and we have 
||u"|| 2 = u' T U' T U'u / . The expanding cone S^, x -^ at x' = T(x) is determined 
by an inequality ||u"|| > 7||u'||, or 

u' t h; u ' > o, (14) 

where H 7 = U' T U' - 7 2 . 

Thus, the required condition DT X (S2) C S^fa;) is formulated in terms of two 
quadratic forms: If the inequality (IT5|) holds, then the inequality (JHJ) must be 
valid too. 

Let us perform a canonical reduction of the quadratic form u /T H 7 u' by a 
coordinate change. Because the matrix U' T U' is real and symmetric, an or- 
thonormal basis of eigenvectors do, di, d 2 , d3 may be chosen. Then, the matrix 
D = {d ,di,d 2 ,d 3 } is a diagonalizer: D T U /T U'D = {A 2 <%}, i,j = 0,1,2,3. 
The eigenvalues on the diagonal A 2 are supposed to be arranged in the de- 
creasing order. In our case there will be one stretching and three contracting 
directions, so, Ajj > 1, Af 2 3 < 1. Let 7 be selected in such a way that A 2 , > 7 2 , 
A 2 2 3 < 7 2 . Under the transformation D the matrix H 7 also becomes a diag- 
onal: 

D T H;D = D T (U' T U' - 7 2 )D = {(A 2 - 7 2 )^}; 
here one diagonal element is positive and others are negativeQ By additional 
dilatation (compression) S = {s^ 1 ^-}, s = \J A% - 7 2 , si )2 , 3 = \Jl 2 - Af^, 
we get 

h; = s t d t (u' t u'- 7 2 )ds = ^, 

cr = 1, 01,2,3 — — 1- The same transformation applied to the matrix H 7 = 
U-^U" 1 -7- 2 , yields 

H 7 = S T D T (U" 1 ' T U- 1 - 7" 2 )DS = {hij}, 

where = hji. 

1 This property is naturally checked in the course of computations at each analyzed 
point of the absorbing domain: its violation would entail an incorrect operation of 
taking a square root of a negative number. The inequalities for eigenvalues of the 
matrix UjU x ensure fulfilment of the condition that a sum of subsets of the linear 
vector space (that is a set of all possible linear combinations of vectors from the 
expanding and contracting cones) is the full 4D vector space: S2 + C2 = V x . 
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A condition iHM for vector c = {1, ci, c 2 , C3} to belong to the expanding cone 
now reads as c T H^c > 0, or 

c\ + c\ + C3 < 1. 

In the 3D space {ci,c 2 ,c 3 } this corresponds to the interior of the unit ball. 
A condition (fT3"l) that preimage of the vector c = {1, ci, c 2 , c 3 } belongs to the 
expanding cone S% is c T H 7 c < or 

3 3 
hoo + Y (^oqCq + h a0 c a ) + Y hapCaCp < 0. 

o=l a,/3=l 

In the space {ci, c 2 , C3} this relation determines the interior of a certain ellip- 
soid. 

The inclusion DT X (S2) C S^fa;) w ^ n ^ e fulfilled, if the ellipsoid is placed 
inside the unit ball. Let us formulate an inequality sufficient to ensure such a 
disposition. We can evaluate coordinates for the center of the ellipsoid from a 
set of linear algebraic equations 

3 

Y^hapcp = -hoo, a = 1,2, 3, (15) 

8=1 

and then estimate a distance of this point from the center of the unit ball: 

P = \/ci + c 2 2 + cl 

After transfer of the origin to the point {ci, c 2 , C3}, the equation for the surface 
of the ellipsoid becomes J2a 8=1 hapCaCp = R 2 , where c a = c a — c a and R 2 = 

Now consider a symmetric 3x3 matrix h = {h a p\. In the diagonal represen- 
tation obtained with an orthogonal coordinates transformation (ci,c 2 ,c 3 ) — > 
(£1, £2, £3) the equation for the ellipsoid surface takes a form li£f + h^i + — 
R 2 , where h, Z 2 , Z3 are eigenvalues of the matrix h. The largest principal semi- 
axis of the ellipsoid is expressed via the smallest eigenvalue: r max = R/ V^mm- 
Now, an obvious sufficient condition for the ellipsoid to be positioned inside 
the unit ball is given by the inequality 

/ = r max + p<l. (16) 

It completes the procedure of the verification of the expanding cones inclusion 
for the point x. 

It may be shown that the above procedure, applied to the points of the ab- 
sorbing domain U with 7 < 1, is equivalent to the verification of the condi- 
tion for the contracting cones in the domain x G T 2 (U) with the parameter 
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i = I/7 > 1: DT-^CVt) c C%? Hx) . It is so because the cones S 7 and C 1 ^ 
are complimentary sets: S* 7 U (7 1 / 7 = V. (Here S 1 corresponds to the cone 
of vectors that either expand or contract but no stronger than by factor 7.) 
Hence, fulfillment of the inequality ([TBI) checked inside U for two parameters 7 
and I/7 would imply that both conditions for expanding and for contracting 
cones are valid in the domain T 2 (U) which contains the attractor Q This is 
sufficient to draw a conclusion on the hyperbolic nature of the attractor. 



5 Numerical verification of the hyperbolicity 

Numerical tests that we performed according to the above method for the 
Poincare map of the system confirms that the theorem is fulfilled. Thus, 
the chaotic attractor under our consideration is indeed hyperbolic. 

For our computer programs we employed algorithms for computation of eigen- 
values of real symmetric matrix, for matrix inversion and for solution of a set 
of linear algebraic equations found in [19]. 

Necessary conditions, required for the expanding and contracting cones to be 
defined, are the inequalities for eigenvalues of the matrix UjU^: A 2 , > 1 and 
A 2 < 1, respectively. Performed computations demonstrate that these inequal- 
ities are valid in the entire absorbing domain for the parameters set ( I12p . The 
the toroidal absorbing domain is defined by Eqs. (flT!) where x an d ip are kept 
constant while and 9 are varied from to 2tt with the step 27r/50. At each 
point we compute the matrix XJ X and find corresponding eigenvalues Aj. Fig. [6] 
represents Aq and A 2 vs. two angle variables at % = 1 (this corresponds to the 
surface of the absorbing domain) and at ip = 0. Observe that Aq > 1 while 
A 2 < 1 (and two rest eigenvalues A^ 3 are also less then one). Similar graphs 
appears at other values of ip and so the same is for the interior of the absorbing 
domain. 

Fig. [7J shows values of the function / ( IToT) that it takes on the surface of the 
absorbing domain atx = l,V' = 0,7=l.l and at the parameters values (fT2j) . 
We see that the function / is always less then 0.5. The function is smooth and 
varies sufficiently slow, so that there are no peaks which can be candidates for 
violation of the inequality (fT6l) between the checked values of <fi and 9. Thus, 
we can conclude that the inequality ffT6l) is valid. 

One can see from Fig. [7] that / depends essentially on only, while the depen- 
dance on 9 is very weak. The dependance on if) is tested to be also weak. Plots 

2 The cones S" 7 and C 7 have a common border only at 7 = 1, while for 7 > 1 they 
do not intersect, as required by the theorem condition: S2 n Cl = 0. 



13 




Figure 6. The first and second eigenvalues of the matrix UjU x vs. two angle variable 
at the boundary of the absorbing area and at ip = 0. 




Figure 7. Verification of the inequality (|16p on the surface of the absorbing domain 
at if) = 0, 7 = 1.1 and at the parameters values (|12p . Observe, that / < 1. 

calculated at different values of ip looks almost identical to the that shown in 
Fig. [7J Qualitatively similar plots appear also at x < 1 that corresponds to 
the interior of the absorbing domain and the inequality (fT6l) remains always 
valid. 

To verify the validity of ([TrJj) for the whole absorbing domain, we fix certain 
value of x an d find the maximum of / as a function of three angle variables 4>, 
if), 6. Thus, varying x, we obtain a function / max (x) that is shown in Fig. [8](a). 
Notice that the inequality / max < 1 is valid both inside the absorbing domain, 
i.e., at < % < 1, and far beyond. 
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12 3 4 5 0.6 0.8 1 1.2 1.4 1.6 1.8 

a) * b) y 

Figure 8. a) Maximum of f(cj),tjj,9) vs. x at 7 = 1.1. b) Maximum of f(4>,ij;,8,x) 
vs. 7. The parameters values are given by (|12[) . 

Now we need to examine the validity of (|T6|) at different 7. We assume that 
/ is a function of four variables <p, ip, 8 and x an d find its maximum at each 
7 with additional requirement < x < 1. As follows from Fig. El^a), / max 
grows monotonically vs. x- Therefore, when a numerical procedure, that seeks 
a maximum, goes out of the absorbing domain and asks a value of / at x > 1 > 
we return / at 2 — x instead. In fact, we create an artificial maximum at the 
boundary of the absorbing domain, that prevents the overrun the boundary. If 
the seeking procedure can not find a maximum in the interior of the absorbing 
domain, it finds the maximum on the boundary at x — 1- The resulting graph 
/max (7) is shown in Fig. EJ^b). We see that the inequality / max < 1 is definitely 
fulfilled in a wide range of 7 both below and above the point 7 = 1. This 
confirms the expected mutual location of expanding and contracting cones 
in the absorbing domain. Hence, we conclude that the hyperbolicity of the 
attractor of the system is established. This result can be reproduced for a 
wide range of values of the parameters of the system . 



6 Conclusion 



In this paper we introduce amplitude equations for a system of two coupled 
non-autonomous self-oscillators that have uniformly hyperbolic chaotic at- 
tractor. We provide an evidence that the system of amplitude equations also 
demonstrates the hyperbolic chaotic dynamics. Given a certain set of param- 
eters, we find a toroidal absorbing domain that is mapped into itself and 
contains an attractor of a Poincare map for this system. Performed compu- 
tations confirm the validity of sufficient conditions for the hyperbolicity. The 
conditions are formulated in terms of inclusions of expanding and contracting 
cones that are defined in a tangent vector space associated with the points of 
the absorbing domain. 

Because of the universality of amplitude equations, the considered system can 
be a model for dynamics of various physical systems. With this example it is 



15 



possible to construct other models with hyperbolic chaos, exploiting structural 
stability of the hyperbolic attractor. A physical experiment demonstrating at- 
tractor of this type has been performed already on a basis of coupled electronic 
oscillators [12]. In applications, the systems with hyperbolic chaos may be of 
special interest because of their robustness (structural stability). An interest- 
ing and now a substantial direction is constructing chains, lattices, networks 
on a base of elements with hyperbolic chaos [6] . Models of this class may be of 
interest for understanding deep and fundamental questions, like the problem 
of turbulence. 

This work was partially supported by DFG and RFBR (Grant No. 06-02- 
16619). 
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